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Long Range Detection by Star Occultation’ 


“8 Harvey Dubner} 


ntroduction 


This paper briefly presents the “Star Occultation” 
schnique for providing long range detection of objects 
1 Space (such as satellites and other space vehicles). 
‘he basic principle involved is that an object moving in 
pace must eventually pass between the observer and 
ome stars; that is, the body occults the light from these 
tars. Thus, the term “star occultation”. The angular 
Osition and time of these occultations are used to 
etermine accurately orbits of objects in space. 

This technique was originally conceived for ground- 
9-air operation; however, it was ineffective in the day- 
me and quite limited for nighttime usage due to atmos- 
heric effects (such as star twinkling). With the advent 
f the space age, an observation platform can now be 
stalled above the atmosphere so that this technique 
wr long range detection by star occultation becomes 
teresting. 

The capabilities, instrumentation and limitations of 
is technique are discussed. Detection ranges of over 
,000 miles are possible; and in one example, 120,000 
iles is shown to be feasible. 


JELD OF VIEW= I710* x 1710° AVERAGE NUMBER OF STARS 


OCCULTED=20 STARS/ DEGREE 


100 FT OBJECT 
100 MILE DISTANCE 
20 INCH TELESCOPE 


CULTED STARS 


110° 


Fig. 1. Detection by Star Occultation 


Figure | illustrates the basic principle and presents a 
aple example of star occultation; in this case, two 
scultations have taken place in yo degree of motion. 


istronomical Data 


The study of star occultation tec hniques began with 
1e compilation of data relating to the number, density 
ad distribution of stars as a function of their bright- 


'* Revised version of paper presented at the 5th Annual 
eeting of the AAS, Washington, D. C., December 1958. 

f Avion Division, ACF Industries, Incorporated, Paramus, 
Pw Jersey. 
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ness. The Milky Way conveniently serves as the galac- 
tic equator with all galactic latitudes chosen in proper 
relation to this reference. 

The relationship between density and galactic lati- 
tude is shown in Fig. 2, a graph of star density (star/ 
degree”) as a function of telescope aperture and star 
magnitude. The 10, 20, 40 and 90 degree curves repre- 
sent the densities observed at these latitudes. The 20° 
curve represents the average density of the overall 
galaxy. The graph illustrates that star density is 
approximately linear with telescope area in the vicinity 
of 13th magnitude stars. To emphasize the information 
obtained in the graph shown in Fig. 2, several significant 
points on the 20° curve were selected and translated on 
to the pictorial display shown in I'ig. 3. 

The average number of stars occulted per degree of 
motion is a function of the effective telescope diameter, 
the linear breadth of the object (transverse to the line 
of motion) and distance. Examples are tabulated in 
Table 1 from the basic formula, 
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| SQUARE DEGREE — 
5" DIA. TELESCOPE 


AVERAGE NO. OF 5 
STARS / DEGREE 


TELESCOPE 
DIAMETER 
| INCH 
S INCH 85 
10 INCH 320 
20 INCH 1,200 


100 INCH 12,000 


Fig. 3. Star Density 


Occultation/degree 


linear breadth of object X effective 
_ telescope diameter 
distance 


TABLE 1 


Average Number of Occultations 


Lins Beat | Distance | Efegve Telescope] Neclons 
200 miles 20 inch 1 
; 100 inch 10 
10 ft. —-- IF 
2000 miles 20 inch all 
| 100 inch | | 
200 miles 20 inch | 10 
100 inch | 100 
100 ft. | = 
| 2000 miles | 20inch (| 
100 ineh 10 


Instrumentation 

The proposed means for implementing the star- 
occultation technique consists of a telescope, a detector, 
a storage and comparator device, a clock, and a com- 
puter. (The computer would most likely be earthbound. 
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Fic. 4. Instrumentation 


Thus a data link for telemetering the data would be 
necessary.) 

Figure 4 shows a telescope mounted on a platforn 
which is space-stabilized by optically slaving it to twé 
reference stars.* A satellite or space station woule 
serve as a suitable platform. Transferring the view fron 
the telescope to a practical electrical output is accom 
plished by an integrating image tube such as an image 
intensifier orthicon. The electrical output is then storet 
for comparison with the next frame to determine th 
occurrence and location of occultations. Occultatior 
location and time are fed to a computer in which orbit 
and trajectories are calculated. 

Telescope—The telescope collects the light energy 
from each star and focuses it on the detector. Thi 
telescope must have a suitable degree of resolution fo 
the accurate identification of a star. Referring to Fig 
3, it can be seen that for a 20-inch telescope stars hay 
an average spacing of about = of a degree. Th 
telescope resolution would therefore not have to b 
much better than 1 minute of are to determine whiel 
star has been occulted. 


The 20-inch telescope normally used by astronomer 
has a resolution of about 4} second of are which fa 
exceeds this requirement. Consequently the foea 
length of the telescope can be decreased, thereby pel 
mitting a reduction in the size and weight of the overal 
equipment. 

On the basis of available data and equipment, i 
appears that an f/1 concentric telescope would b 
sufficient to achieve the desired resolution over larg 
fields of view. 

Detector—The detector is used to convert the ligh 
energy to its electrical equivalent. Investigations ¢ 
available devices for converting light to its electrice 
analog have indicated that the image intensifier orthi 
con developed by RCA fulfills the essential require 
ments of a detector. Its potential sensitivity is si 
times better than the human eye, and its resolution 1 
100 lines per inch, which is adequate. Present models ¢ 
this orthicon have a sensitive surface that measure 


* Only two stars are required to establish a reference. 
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Fie. 5. Image Intensifier Orthicon 


6 inches x 1.6 inches. However, RCA states that 
ger sizes can be produced. 

Figure 5 illustrates the basic operation of a detector. 
tarting with the light energy striking the photo- 
missive surface (cathode), electrons are emitted which 
re attracted to the intensifier due to its high potential. 
ach electron, in bombarding the intensifier surface at 
igh velocity, dislodges several other electrons from 
1is surface. These secondary electrons, in turn, strike 
1e target surface. Scanning is then performed on the 
ther side of the target surface in the same manner as 
} an ordinary orthicon and amplified in a conventional 
ype of multiplier. Internal noise in this orthicon is low 
aough so that individual photons can be counted. 

For optimum sensitivity, frame time should be 
pout equal to the time of occultation. Considering 
ve case of a 10-foot object traveling at an orbiting 
elocity of 20,000 feet per second, occultation time would 
® spy second. This presents a problem which would 
ave to be solved by the use of several detectors in 
arallel. 

Storage unit and comparator—Determining when an 
scultation has occurred would be accomplished by 
ame-to-frame comparison of star locations within the 
lescopic view. Considering the nature of the necessary 
ita, the use of electronic storage tubes seems to be the 
gical choice. 

Storage tubes capable of handling up to 250,000 com- 
etely isolated bits are now available; writing and 
ading rates up to 500 frames per second have been 
ed. Two storage channels are being considered for 
is system—one for writing the incoming information 
1d the other for comparison with it. Comparison 
ould then be effected by the process of simple sub- 
action, with the resulting occultation information 
sing fed to the computer. 

Computer —A digital computer would be employed 
r the calculation of orbits and trajectories. Since the 
‘tical resolution selected is consistent with star spac- 
g, information fed to the computer will be accurate 

time but only approximate in angular data. Hence 
e computer must store information of star locations 


within its memory if it is to provide angular data with 
an accuracy in the order of fractions of a second of are, 
as are desirable in computing trajectories. 

As an example of the accuracy attainable, consider a 
satellite in a circular orbit as it passes directly over- 
head. Only two occultations are required to derive its 


TABLE 2 
Occultation Technique Accuracy 

h = Distance from earth.......| 200 mi. 2000 mi. 
Time for one revolution......... 91 min. 156 min. 
Maximum time orbit is visible 

from earth’s surface.......... 9 min. 41 min. 
Time between occultations...... 15 see. 26 sec. 
Ah due to time error of .01 sec. 

Of ities eisai Sasi 5.5 mi. 4.6 mi. 
Ah due to angular error of 1 sec. 

OLSATCrerseust CS Pen ee: 2.33 mi. 3.30 Mi. 


Assumptions: 1. Circular orbit around earth; 2. One star 
occulted per degree of are. 


orbit. Table 2 contains the data for two targets at dis- 
tances of 200 and 2000 miles above the surface of the 
earth for the conditions of two occultations separated 
by 1 degree. Input errors of 0.01 second in time and 
1 second of are have been assumed. The orbit errors 
are capable of being reduced by additional occultations. 

The knowledge that a target is following a ballistic 
course may be utilized to improve system sensitivity. 
When operating at low signal-to-noise ratios, a false 
indication of star occultation may occur. However, if 
this occultation were stored and compared with future 
occultations, it could be evaluated as to whether it was 
false. In effect, this is equivalent to a very powerful 
scan-to-scan correlation technique. By increasing the 
computer capacity, operation at lower signal-to-noise 
levels is feasible. 


Limitations 
General 

Limitations of the star occultation technique stem 
directly from the nature of light itself. They fall into 
two categories in accordance with the dual character 
of light: (1) the wave character of ight imposes range 
limitations; (2) the granular character of light imposes 
velocity limitations. 


Range Limitation 


An object passing in front of a star does not cast a 
true geometric shadow—it forms a diffraction pattern. 
The smaller the object or the greater the distance to 
the observer, the larger the diffraction pattern, until 
finally the object cannot be discerned as occulting a 
star. Table 3 contains a chart indicating the ranges and 
object sizes at which this limitation begins to be notice- 
able. 

It should be realized though that the ranges shown 


The Journal of the Astronautical Sciences 3 


TABLE 3 
Diffraction Limitation 
i i R for Reduction of 
LG vans ge Ne 2 8 Effective Occultations 
Feet Inches Miles 

10 20 2,400 

100 12,000 

100 20 24,000 

100 120,000 


R « (Effective Telescope Diameter) 
x (Effective Target Diameter) * 


* Assuming spherical target. 


represent a decrease in detection probabilities and not 
complete absence of information. A bright star will still 
give evidence of a partial occultation. Of interest is the 
fact that available equipment will effectively operate 
right to this limit. Astronomers obtain data by working 
into the diffraction pattern. 


Velocity Limitation 


General—Velocity limitation, the more serious of the 
two limitations, is due to the granular nature of light. 
Light from a star travels in small packets of energy 
called photons. The ability to detect a star is deter- 
mined by the presence of a certain number of photons 
per second striking the detecting device. In the illus- 
trations used thus far, 16,000 photons per second are 
entering the telescope aperture for the dimmest detect- 
able star. Thus, if an object occults a star for about 
Teo00 Second, it would be difficult to determine whether 
the star has been occulted, since in such an interval the 
chances are about 37 percent that a photon would not 
have arrived. 

In addition, the characteristics of the detector create 
further difficulties since the detector does not neces- 
sarily respond to each photon. The detector is only 
responsive to certain wavelengths, and then only to a 
percentage of photons within that wavelength. For 
example, the efficiency of the eye to sunlight is about 
1 percent, which means that, on the average, 100 
photons are required to obtain a response from the eye. 
The image intensifier orthicon has an efficiency of ap- 
proximately 6 percent. 

False Alarm Rate-—Due to the fluctuating photon 
count, the problem exists of determining whether an 
occultation or a “false alarm” has occurred. In Appen- 
dix I is an outline of the method for determining the 
average false alarm rate. 

Figure 6 contains a plot of average false alarm rate 
vs. occultation time for 16,000 effective photons per 
second from the dimmest detectable star. This figure 
indicates a rapid decrease in the false alarm rate with a 
very small increase in occultation time. Therefore in 
order to obtain a reasonable false alarm rate, the time 
of occultation is generally such that an average of 12 
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100 


FALSE ALARMS PER HOUR 


0.1 
72 4 6 8 


OCCULTATION TIME (MILLISECONDS) 


4 8 12 16 
AVERAGE PHOTONS PER OCCULTATION 


Fic. 6. Average False Alarm Rate vs. Occultation 
for 16,000 Effective Photons per sec. 


t 
to 16 photons normally arrive from the dimmy 
allowable star. ' 

In the preceding illustrations indicating the numb 
of occultations per degree with various telescope a 
object sizes, allowances were not made for veloei 
limitation. Table 4, however, contains a tabulation 
the velocity at which effective occultation decreas 
assuming an average false alarm rate of 1.0 occultati 
per hour per star. 


TABLE 4 
Velocity Limitation 
: Photo- Average Obj oe 
Detector Electron Effective — 5 we a 
E Photons/Sec Velocities a 
| Feet Miles per hon 
Ideal | 1.0 16,000 10 7,000 
100 70,000 
Image Intensi- | .06 1,000 10 550 
~ . | 
fier Orthicon | ) 100 5,500 
Eye (Equiv.) eae (i 160 | 10 100 
100 1,000 


Detection probability —As a target sweeps across t 
observer’s field of view, normally three occultatio 
will provide a reasonably accurate determination of 
path. However, because of the erroneous informati 
produced by false alarms, four occultations are deem 
a practical minimum for obtaining reliable and accure 
data. Hence the problem becomes one of determini 
the star density required to ensure a high probabili 
of obtaining at least four occultations during the ti 
that the target is within the field of view of the te 
scope. Figure 7 graphically depicts this probability 
a function of average number of occultations. T 
average number of occultations required is ten for 
percent detection probability. 
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Fig. 7. Probability of Occulting 4 or more Stars 


Whether a star is usable depends on the time of 
ccultation which, in turn, is dependent on the length 
md velocity of the target. These factors are shown in 
heir relation to telescope diameter in the following 
quations: 


1 
Occultation Time 


‘elescope Diameter « ( 


(1) 


1 pe 
es Angular Width of Tesi) 
r, 


Velocity 
Length 


; 1/2 
y, eae) (2) 


Teles 1 ox 
elescope Diameter ( Width 


These equations coupled with a knowledge of star 
ensity and energy enable one to calculate the size of 
ie telescope required for a given set of target condi- 
ons. A typical case is shown below: 


Target = 20 feet x 20 feet 
Range = 200 miles 
Velocity = 18,000 mph (Orbit near 
earth surface) 

False Alarm Rate = 100/sec 
Detection Proba- = 99% 

bility 
Telescope Diameter = 31 inches 
Field of View = 40° x 40° 
Number of Effec- = 360,000 

tive Stars 


bility to Detect Orbiting Vehicles 


To illustrate the potentiality of the star occultation 
chnique, consider a vehicle in orbit about the earth. 
quations 3, 4, and 6 contain a few basic relationships 
hich emphasize the desirable fact that range varies as 
ie fourth power of the telescope diameter and the 
urth power of the target diameter. Equation 3 is 
mply a restatement of an equation previously de- 
loped; equation 4 is a statement of the relationship 


that exists between target velocity and distance to the 
center of the earth; equation 5 is a simplifying assump- 
tion; and equation 6 is the result derived from these 
three parts. 


(Velocity )”(Range)"” (3) 
(Target Diam) 


(Telescope Diam) « 


For Target in a Circular Orbit 


: (4) 


Velocit 
(Velocity) « (Distance to Earth Center)! 


For Simplicity, Assume 
Distance to Earth Center = Range (5) 


Thus, for Photon Limitation 
(Telescope Diam)4(Target Diam)! « (Range) (6) 
Detection capability is demonstrated in the tabula- 


tions below, in Table 5, which are the results obtained 
with Eq 6. Also included is the diffraction limited range. 


TABLE 5 
Detection Capability 


Photon Limited ; . i 
Target Pelescons kanes for Dee es 
Feet ‘ Inches Miles Miles 

100 30 5,000 34,000 
40 16,000 48,000 
60 80,000 68,000 
30 100 5,000 34,000 
40 16,000 48,000 
100 620,000 120,000 


* For ease of calculations: 1. Target assumed to be in circu- 
lar orbit; 2. Distance from target to earth center assumed 
equal to range. 


APPENDIX I 
Derwation of the False-alarm Rate 


f = false-alarm rate for a particular star 
F = average false-alarm rate per star for all stars 
N = number of effective photons arriving at the 
telescope aperture in | second from the star 
under consideration 
T = the occultation duration 
At = the trial time 
p = the probability that no photons arrive in the 
trial time At 
q = the probability that a photon arrives in the 
trial time At 
r = T/At, the number of trials during an entire 
occultation 
uw = mean number of trials required to produce a 
run of no photons of length r 
S(N) = total number of stars which produce a photon 
intensity of at least N at the telescope 
aperture. 
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From standard probability theory, * 
= ieaiben (haat (7) 


qp" eas) 


U 


Assuming a Poisson distribution for the arrival of 
photons from a particular star, the probability that no 
photons arrive during At is 


p= Bet (8) 


Substituting this distribution and the definitive form 
of r into Eq. 7 one obtains 


TN 
Apres aaeae (9) 


1— ¢e%4t 


For small At, Eq. 9 may be rewritten as 


TN 
e —1 


NAt 


The false-alarm rate per star is the reciprocal of the 
mean recurrence time for a run of length r: 
1 


f= ete Nike se (11) 


uU= (10) 


The average false-alarm rate is the summation of the 
false-alarm rates of all the stars divided by the total 
number of stars under consideration. 

The summation extends from the lowest-intensity 
star of interest to the highest intensity stars, which for 
practical purposes can be considered of infinite in- 
tensity: 


ee (12) 


where AS is the number of stars which have the false- 
alarm rate f. In the limit, 


ioe) 


i| 
(if = en ? (13) 


~ S(No) 


* W. Feller, Probability Theory and Its Applications (John 
Wiley & Sons, Inc., New York, N. Y.) p. 286. 


where S(No) is the number of stars which have intensi 
exceeding a given intensity. 
From Fig. 2, a relationship can be derived for S( 
as a function of the photon intensity: 
3 
= — 1 
S(N) N 

The relationship is reasonably accurate for star ma 
tudes brighter than +15, the region of most intere 
Differentiating and substituting in Eq. 13, 


Letz = TN; then 


F = —No 


TNo 


For reasonable false-alarm rates, T’No, the averag 
number of effective photons blocked during an occulta 


tion, is rarely less than 10. Therefore, 
eo i (17 
i 


so that Eq. 16 can be closely approximated by 
ioe] Ca “ 
F = —No if ae (18 

TNo v 


This can be expanded in an asymptotic series giving 


IN 
e 0 


SS : ] 
F 7 (1g 


Equation 19 is plotted in Figure 6. 

It is interesting to compare the average false-alart 
rate, Eq. 19, with the false-alarm rate of the dimme: 
star under consideration, Eq. 11, 


FCN) eel 
f(Nc)  TNo 


(20 


TN,» is about 12 to 16 in most of the examples used. 


ANNOUNCEMENT 


The Journal of the Astronautical Sciences welcomes papers on any aspect of astronautics. 
Contributions should be original, generally quantitative in nature, and should satisfy high 


standards of scholarly excellence. Papers are especially solicited in the fields of space flight 
mechanics, space vehicle design, space physics, propulsion, guidance and control, communi- 
cation, space medicine and astrobiology, and applications of astronautical systems. How- 
ever, any other papers concerned with astronautical investigations will be considered. 
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~ On the Problems of Re-entry into the 
Earth’s Atmosphere’ 


Alfred C. Robinson and Algimantas J. Besonis 


bstract 


Re-entry into the earth’s atmosphere has been studied 
om the standpoints of deceleration, heating, and accuracy 
impact. This has been done for re-entry speeds consistent 
ith return from near satellite orbits, and for speeds con- 
stent with re-entry from a circumlunar orbit under several 
figurations of lift and constant or variable drag coefficient 
sumptions. Heating considerations are based only on stag- 
ition point influences. It is shown that deceleration and 
ak heating rates are not larger than those occurring in 
ullistic missile re-entries. The total heat input, however, is 
uch larger as the heating occupies a much longer time. It 
ypears that simple, non-lifting re-entry will be feasible from 
tellite orbits. The lunar re-entry, on the other hand, pre- 
nts a severe total heat problem and accuracy requirements 
e such that some lift or other control will probably be 
quired. 


EC I—Introduction 


Before studying the re-entry problem, it is well to 
ate the criteria of a successful re-entry. The mission 
ing considered here is that of recovering an instru- 
ented or manned vehicle after it has been orbiting 
e earth or gone around the moon. There are two parts 
this problem: it must reach the surface of the Earth 
ithout substantial damage to crew or payload; and 
must be found. The first aspect is concerned with 
rodynamic heating, deceleration, and dynamic pres- 
re; the second with impact point prediction. 

Very little will be said about the configuration of the 
hicles considered. This is a separate field of investiga- 
iy and an important one. It is assumed that the bodies 
e generally blunt and are hemispherical in the region 
the stagnation point. It is further assumed that the 
xgnation point heating is the most severe on the body, 
d that therefore it is the only point for which the 
ating needs to be calculated. 

Three types of re-entry are considered: (1) no lift, 
nstant drag coefficient, (2) no lift, variable drag 
efficient, (3) low lift. All the studies reported herein 
» concerned with the behavior of the vehicle inside 
e atmosphere. All simulations and computations 
vered deal with motion, heating, etc., from ‘“‘re-entry 
itude” down to the surface of the earth. Above the 
entry altitude, the missile motion may be computed 
sed on vacuum conditions, except in certain special 
ses, such as lifetime computations. 


* Aeronautical Research Laboratory, ARDC, August 


8. 


Fia. 1. Coordinate System 


In an exploratory study such as this, high accuracy 
is not required per se. Therefore it was decided to use 
an analog computer and take advantage of the high 
speed and flexibility of that equipment. This made 
necessary one or two points of procedure which could 
have been omitted otherwise, but it is felt that the 
utility of the analog computer was such as to justify 
some slight additional analytical complexity. 


SEC If—Development of the Equations 


In this section the equations are developed which 
were used to describe the system mathematically. 
They fall logically into two divisions: mechanical equa- 
tions of motion, and the heating equation. They make 
use of different techniques and disciplines and are 
largely independent. 


A. The Equations of Motion 


It appears to be adequate for present purposes to 
consider motion only in a single plane. It would appear 
that all the important effects can be studied except 
possibly impact errors. These can be treated as per- 
turbations from the plane. Furthermore, the equations 
will be developed for a non-rotating earth. 

Figure 1 shows the definition of the quantities to be 
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studied. The equations describing this motion are 
given below: 


v= a(S) —o(1- Ben 


2h SC, 
y= -£(1 2) cosy 6+ 4 (5) 


bh = V sin 
e (1) 
_ 24 _ 
b=E(1- Py eos. 
r=Reth, 


It will be observed that a first order power series ap- 
proximation has been applied for the terms involving 
the reciprocal of r and the reciprocal of r?. This will 
result in an error less than 0.4 per cent in computing 
the gravity force. The dynamic pressure is being com- 
puted from a differential equation rather than explicitly. 
This technique has been set forth by the senior author 
in a previous work. [4] 

The quantity K comes from the assumed exponential 
variation of the atmosphere. It has been assumed that 
the density may be expressed by 


p(h) = pye*" (2) 


where po and K are adjusted to give a best fit for the 
actual atmosphere desired. The values which were 
chosen were 


po = 2.78 X 10-* slugs/ft*, and 
K = 4.276 X 107° ft7. 


This fits the ARDC model atmosphere to about 20 per 
cent from sea level to 400,000 feet, which is the altitude 
range of interest. It appears that this is accurate enough 
for studies of this kind. 


B. Aerodynamic Heating 


Only the region of maximum heating (the stagnation 
point) is considered since it represents the most critical 
design area. Blunted bodies of radius R will be con- 
sidered, since only these types of configurations appear 
to be practical, because of the high aerodynamic heat- 
ing rates encountered under re-entry conditions. Since 
the maximum aerodynamic heating for both re-entry 
conditions occurs below 300,000 feet altitude, the exist- 
ence of chemical equilibrium in the stagnation region 
will be assumed. 

Experimental results of Rose and Stark [5], show 
good agreement with heat transfer rates at the stagna- 
tion point of blunt bodies at hypersonic speeds assum- 
ing an equilibrium boundary layer, as stated above. 

Under these conditions the stagnation point heating 
rate is given by Fay and Riddell [6). 
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Qa = 94g (puhtw)°"\(Pekts) {1 = (Lz = 1) (ho/hs)} 
{(hs — hw) (du/dx)s”} 
assuming a Prandl number of 0.71. 

Equation 3 in its present form is rather lengthy 
analog simulation. Order of magnitude analysis i 
cates that it can be simplified and reduced to a relat 
lending itself more readily for the present ‘pu 
while still retaining the essential parameters influ 


ing stagnation region heat transfer..__ 
Consider the following approximations 


Clee ae r ee 1 


(pu )we ~ (ou)! 


Under these assumptions Eq 3 reduces to 


2 
Qa = .94g( pu)?” S (du/dx)s° . 


et 4 FZ Oh eee, 


An approximate solution for the nondimensiot 
velocity gradient at the stagnation point in hypeam 
flow is given by Ting Yi Li [7] as 


rac Mars\cew.7) 


Then assuming that a Sutherland type viscosity + 
efficient relation holds 


=9 po F(k) 
Q. = 6.2076 X 10 E a Ue ae 


ee ee 


F(k) is shown in [9] as a function of k. 

The density ratio k has been computed by Feldm 
[8] as a function of velocity and altitude up to satell 
velocities for chemical equilibrium behind the she 
wave. It should be noted that the density ratio k at - 
high velocities considered equals the density ra 
behind the shock wave consistent with the assumpti 
made by Ting. This equality is well indicated 
Feldman. 

The density ratio k for lunar re-entry bodies (ini 
velocity of 36,000 ft/sec) can be computed.from tak 
computed by Gilmore [9], for equilibrium composit 
of air.and the normal shock relations in a man 
similar to the computation of Feldman. 

Computations of the parameter [F(k)/k]®> for 
configurations considered indicate that the variat 
of this parameter in the range of high heating rate 
sufficiently low to justify its considerations as a ¢ 
stant equal to 3.7; then Equation 9 reduces to 


0.5 
O5=.23 < 10m: || cae ( 


It should be noted that the assumption of chem 
equilibrium for all re-entry conditions will give sligl 


wer heating rates for the high drag configurations, 
it the results should still be within the accuracy of 
| per cent attainable in experimental verification. 
The assumption of viscosity variation with tempera- 
re according to a Sutherland type relation will give 
wer viscosity coefficients than predicted by Hansen 
a for chemical equilibrium up to temperatures of 
00°K. At higher temperatures a very sudden drop in 
e ratio of the viscosity coefficient as predicted by 
ansen to the coefficient considered is observed, but 
is effect of the variation of the viscosity coefficient 
ith temperature on the heating rate is weakened by 
e square root dependence of the heating rate on the 
scosity coefficient. 


_ Approximate Solution 


As an introduction to the method which was used to 
cast the equations for computer use, it is instructive 
consider an approximate analytical solution which 
is presented by Allen and Eggers [1]. The develop- 
ent which follows is not in any fundamental way 
fferent from theirs. It differs only in the way solutions 
e presented and in the use of non-dimensionalizing 
ctors. 
If the missile enters with sufficient steepness into the 
mosphere, the decelerations will be very large, and 
i much distance around the earth will be traversed. 
nder these conditions, the earth might be considered 
flat and non-rotating. Equating force along the rela- 
re wind axes gives for Eq. 1 


V = —po et oe oe + gsin y 
(11) 
y= af COS ¥. 
is also necessary to use the relation h = V sin y to 


mplete the set of equations. It is assumed that 
vity is a negligible force and that SCp/m is constant. 
th these assumptions are relatively good at the high 
tudes and high speeds which are usually of interest. 
gravity is negligible, the flight path angle will not 
ange. Negative flight path angles will be of interest 


we substitute y = —6 and the equation becomes 
as = —Vsin@é 
dt 
(12) 
dV Po SCp —Kh 72 
———— sil 
dt 2m 


iminating time between these two equations gives 


dV nee SC 


eae ah. (13) 
V 2m sin 0 


this point we make a change of variable 


aK} = ln ey ey (14) 
m 


2K sin @ 


With this change, Eq 13 becomes 


a = 6 dz (15) 
which may readily be integrated to give 
V —e 2 
re = fle). (16) 


If Eq 14 together with Eq 16 are substituted into 
the expression for acceleration in Eq 12, it may be 
shown that 


d : ee 
= = —Ksin@ Veere-.. (17) 
It may easily be shown that the function e* ¢” ~ 
has a maximum value of 1/2e and that this maximum 
occurs at « = 0.693147. Therefore, the maximum value 


of acceleration is given by 


dV K sin 0V x" 
fe te 


so the ratio of the acceleration to the maximum accelera- 
tion is 

. dV 

dt 


a) z 
dt max 


This is the second of the functions to be tabulated. 
From Eq 12 it may be seen that 


dh dx 


Deh te oleh aa) (19) 


LY es Oe 2 
ai V sin @ = Kd: (20) 
From this, it follows that 
e de = — KVzsin 6 di. (21) 


In integrating this expression, a question arises as to 
the definition of zero time. Arbitrarily, the time corre- 
sponding to maximum deceleration has been selected 


. . . . 2 x 
as the origin of time. The integrand e° dx may be 
transformed by substitutions of e~* = z. Using this 
transformation 


fru -[ 


ae ar a (22) 


at the point of maximum deceleration, e~* = }. There- 
fore, the integral to be evaluated is 
a aE) 
e ae. & 
/ 2 Vee Si Obs re (23) 
05 2 


where 7 is a non-dimensionalized time. Evidently, 
Eq 23 may be written 
rene 0.5 


i ie ede ike a = felt.) (24) 


Zz 


This is the third of the tabulated functions. This in- 
tegral has been tabulated by Blanch. [11] 
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There is another set of equations which is important 
in this problem, namely those relating to aerodynamic 
heating. It has been shown by various authors that the 
aerodynamic heating rate at the stagnation point may 
be expressed by an equation of the form 


yee wear (25) 


If V is in feet per second, p. in slugs per cubic foot and 
R in feet and Q, in BTU per square foot per second, 
and a = 3, then A is approximately 2.3 times 10~*. 
It is possible to non-dimensionalize Eq 25 in much the 
same way as the earlier equations and the result is 


2K sin 0 


5 A ay aa ae —x/2 —ae—* 
Qu a / CDS Vie ee (26) 
VR = 
m 
It may easily be shown that this has a maximum when 
« = In2a and that this maximum is 
J 2K sin 6 


Qaln = ae (5) Lica (27) 
m 


so that the ratio of the heating rate to the maximum 
heating rate is given by 
Qu 2. p~%/2—ae7 
— = /2a¢ e ae een) (28) 
One 
This is the fourth of the tabulated functions. Equation 
26 may be integrated over x to give the total heat input 
down to any specific value of x. The result is 
ieee ST 2 


VR K sin 0 C 2 
m (29) 


zx 
—«/2—(a—1)e-* 
| ee ee 
=-s0 


This integral may be transformed by use of the sub- 


stitution (a — 1)e~* = w2, so that Eq 29 becomes 
AVa 2 
Q; = ~~ = may) ae 
VR _ Cs 
K sin 6 
m (30) 
9 Vanie 2/2 . 
tls = Caan du 
a/ a— 1 0 
The final value of Q, is 
AVS” L 


a 


Qal i 


ss ; [or 
JR Pate on 
KO sin @ (eee 
m 
The ratio of the total heat input to the final value is 
given by 
Qa 2 ie cee 
= e du = f(a, 2) 32 
Ol, veo A tales. 


which is the fifth of the tabulated functions. This 
integral is simply the probability integral and has been 
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tabulated by numerous authors. To summarize the 
tabulated functions are listed below: 


fla) == 
fo(a) = _— Dota 
et of 
fs(a) -r=[ © dz 
4) 4(d, x) = if 23 /2ae ota 
Qa ola 2 anne AM 
fs(a,z) = AEA a Fae ve 


fi(x), fo(w) and fi(x) were computed to six-fig 
accuracy using tables given by Comrie [12] and roun 
off to four significant figures. f;(x) was tabulated us! 
the nine-figure tables of Blanch. The probabil 
integral fs(x) was taken from Peirce [13] though 
other standard table of the probability integral wo 
serve. 


TABLE 1 : 
Values of Functions * 
x fila) Sox) F3(x) f4(3, x) RCH 
. 
0. 0. 0. 0. g 
S50) 0660 0643 0019 9990 
58 1080 1412 0076 9971. 
—0.6 1617 2590 +3.8703 | 0230 9931 
—0.4 2250 4104 +2.8226 | 0562 9854 
=i 2948 5771 +2.0472 1144 9729 
0.0 3679 7368 +1.4409 2011 9545 
0:2 4410 8656 +0.9451 3134 9297 | 
+0.4 5115 9536 +0.5247 4426 8985 
+0.6 5776 9955 +0.1573 5766 8617 
+08 6381 9945 =f. TTAS 7032 8200 
+1.0 6922 9585 —0.4722 8124 7749, 
ai 7399 8965 =. 7515" | Bor 7276 
se es 7815 8187 —1.0143 9570 6794 
+1.6 8172 7330 —1.2644 9902 6312. 
HettS 8476 | 6457 —1.5046 9999 5839 | 
2.0 8734 oo)» BOIS At cee eeners 9899 5381 
ups 8951 4827 —1.9631 9641 4944 
ed ams | 4114 —2.1844 | 9266 4531 
2.6 9284 3481 —2.4009 8808 4148 
2.8 | 9410 | 20927 | —2.6154 8298 3781 
aead0) 9514 2450 —2.8266 8109 3446 
58)0 9601 2043 —3.0358 7215 3136 
tea 9672 1697 —3.2434 | 6675 2852 
+3.6 9730 1406 —3.4495 6150 2591 
+3.8 | 9779 1163 —3.6545 5648 2352 
+4.0 | 9829 0962 —3.8585 5189 2134 
4 9851 0791 —4.0619 4728 1935 
+4.4 | 9878 0651 —4.2647 4313 1754, 
+4.6 9900 0536 —4.4668 | 3929 1589 
+4.8 9918 0440 —4.6687 3574 1440 
+5.0 9933 0360 Se YAY. 3249 1304 
45,2 9945 0297 =5.0r is 2950 1181 
+5.4 | 9955 0243 —5,2734 2678 1069 
+5.6 9963 0200 —5.5038 2429 0968 
+5.8 9970 0164 —5.6769 2202 0876 
+6.0 | 9975 0134 —5.8771 1996 0793 


In order to use this table, it is necessary first to com- 
ute the maximum value of each of the functions. 

Once these have been computed, it is only necessary 
0 multiply them by the appropriate value of function 
rom Table 1 in order to get the value of the variable 
or the associated values of x. It remains, then, to 
letermine the relationship between x and h. This is 
siven in Eq 14. By this means a complete solution 
subject to the above named approximations may be 
ound for any specific re-entry. 


). Normalized Equations 


Consider the following changes of variable: 


K sin V2 
V=V a eet clase Ae 
aad KVssnd * (“ jp 
2e 
m 
| 2K sind Ve 
Q, = 0. —= Fite Siac pene (33 
V/R CS /6e ) 
m 
CyS Po 
eK ire : 
teegst i Ce 


f these transformations are introduced into the Eq 1 
ind normalized Eq 26, the result is 


ee ereh\.. 
ee 55 (: re) sin 


d —B(1— 7A) os7 4 7 


dr Rez dr 2ev \Cd 
dat =C'v sin y Xo = 8.601 (34) 
dr 


dg 2q¢dv _dz do ( h 
ON eee Les ‘pee 5 
dr ovdr q dr’ dr He Rr pet 


1 2 an 
h= = ie + FP} O.= V8V/ar 
yhere 
Bye tig pele Robawey ls ty 
KV? sin 6 So Say se K sn 6 Re 


ee CyS Po 
es ( m ) 2K sin 6- 


‘hese are the equations solved on the computer. 
yhanges of variable were based on the maximum values 
ccurring in the simplified solution of the preceding 
section. The maxima which actually occur during re- 
atry may differ from those calculated. It should also 
e understood that the parameters Vz and 6 used in the 
ormalization are not necessarily the actual entry speed 
nd angle. It is convenient to use values of V» and 6 
hich are near the actual values. 

The “re-entry altitude” at which “re-entry angle”’ 
defined is that corresponding to x = 8.601. From 
‘able 1, it is seen that, at this point in the approximate 
lution, about 0.02 per cent of the velocity loss due to 


TRANSFORMED TIME T 
= E, 


w v ° vt z +4 
NON- DIMENSIONAL ALTITUDE x 


Fig. 2. Normalized Variables vs Altitude 


drag has occurred, the heating rate is about 5 per cent 
of maximum, and about 2 per cent of total heat input 
has occurred. This was selected as the starting point for 
all computer runs though it does not correspond to the 
same altitude in every case. The initial varies from 
300,000 feet to 480,000 feet, depending on re-entry 
speed and angle. The exact value in any particular case 
may be computed from Equation 14. 


SEC III—Re-entry from Satellite Orbits. No 
Lift, Constant Drag Coefficient 


Many different combinations of parameters were 
studied, and while it might be possible to keep record- 
ings such as Fig. 2 for each case, the resulting volume 
of data would make intelligent comprehension 1m- 
possible. Accordingly, for the most part attention was 
fixed on maxima of the various curves. Figure 3 shows 
plots of these maxima for the cases studied. Three 
different re-entry velocities were considered: 25,000, 
26,000 and 27,000 ft/sec. Figure 3A shows the peak 
decelerations, as a function of re-entry angle. Simu- 
lator results for the three re-entry speeds are indicated 
by the solid curves. The dashed curves represent the 
values predicted by the approximation of section IIC. 
It may be seen that in all cases the simulator solution 
approaches the analytical solution as the angle 
increases. The analytical solutions are all straight lines 
passing through the origin. At small re-entry angles, 
the simulator results deviate seriously from the ana- 
lytical results, as would be expected. However at an 
angle of 10°, the analytical result gives a fairly good 
representation. It may be seen that at small re-entry 
angles, the lower the speed, the higher the deceleration. 
At higher re-entry angles, this trend is reversed. 

It is also necessary to know the range over the earth 
which the missile traverses between re-entry and 
impact. Equally important is the sensitivity in range 
to errors in the re-entry conditions. Figure 4 shows the 
range as a function of re-entry angle for each of the 
three speeds studied. The range increases rapidly for 
small re-entry angle. The curves shown are independent 


The Journal of the Astronautical Sciences 11 


PEAK DECELERATION (q's) 


6 C) Ww 
RE-ENTRY ANGLE (DEGREES) 


Fig. 3a. Peak Deceleration vs Re-en- 
try Angle. 
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Fic. 3b. Heating Rate Parameter vs 


Re-entry Angle. 
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of CpA/m , but it should be recalled that this is the 
range from re-entry until impact, and that re-entry 
altitude depends strongly on CpA/m. In determining 
the impact point, this should be kept in mind, 

Possibly more important is the sensitivity of this 
range to various errors. Some of these are shown in 
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Fia. 3c. Total Heat Parameter vs 
entry Angle. 


6 8 Lo) 2 


05 3 1000 

ne & 
el: 
$0 3 600 
z « 
~ a 
3 : 
g 0.3 = 600 
z $ 
= u 
= “ 
Cees 3 W400! 
a = 
% ou x | 

ver PS 

— 
2 4 6 6 10 RE-ENTRY ANGLE (DEG) 
Fie. 5. Range Sensitivity to Re-entry Speed + Angl 
Error. 


Fig. 5. The sensitivity to re-entry angle error is showi 
for each of the three speeds studied, and the sensitiv 
ity to re-entry speed error is shown for 26,000 ft/se 
only (dashed curve). Observe that at about 3°, all th 
curves start increasing rapidly for smaller angles. I 
possible, the re-entry angle should be kept larger thai 
3°, to improve impact accuracy. 

These sensitivities might be made more meaningfu 
by taking an example. It is assumed that the re-entry 
angle uncertainty will be 0.1°, the re-entry velocity 
uncertainty 20 ft/sec, and that the density uncertainty 
is 30 per cent at all altitudes of interest. If thes 
assumptions are used for the errors, the results 0 
Fig. 6 are obtained. The three error contributions ar 
shown as a function of re-entry angle, and also th 
combined error. Again it may be seen that from th 
accuracy standpoint, re-entry angles smaller thai 
about 38° are very undesirable. At small re-entr 
angles, the angle uncertainty is the predominant erro! 
At steeper angles, the density uncertainty become 
dominant. The velocity error contribution is negligibl 
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Fie. 6. Miss Distance for Satellite Re-entry 


; all angles. The nominal re-entry speed assumed for 
ig. 6 was 26,000 ft/sec. ’ 

All the foregoing data are summarized on Fig. 7. For 
re-entry vehicle of this type, there are only a few 
voices open to the designer. The re-entry speed is 
xed by the type of orbit from which the vehicle is 
sscending. There remain only the re-entry angle, the 
rag-mass parameter and the stagnation point radius 
curvature to choose. In Fig. 7, the re-entry angle is 
otted against the product of R and CpA/m. Every 
vecific set of choices of the three variable parameters, 
en, will correspond to a point in the plane of Fig. 7. 
hile re-entry angle is the primary ordinate, peak 
seceleration and impact uncertainty are plotted as 
Iditional ordinates inasmuch as they are single-valued 
nctions of re-entry angle. In this, the re-entry speed 
s been assumed to be 26,000 ft/sec, which is close 
the re-entry speed for the lowest satellite orbits. 
We are now in a position to discuss the problem of 
-entering the earth’s atmosphere with this type of 
hicle. The heating problem is the most complicated, 
it will be considered first. Obviously the most desir- 
le heating situation is one where all the heat which 
ters the body may be re-radiated. This would require 
at peak heating rates be held below about 20 BTU/ 
/sec, and this in turn would require that R CpA/m 
on the order of 200 even for small re-entry angles. 
a 5 ft radius of curvature is used, CpA/m would 
ve to be 40, which means a very light, high-drag 
ucture indeed, less than 1 pound per square foot of 
uivalent drag area. As was mentioned earlier, the 
al heat input curve for this case is of little interest, 
eause all the heat is re-radiated. 

A more typical blunt configuration would have a 
A/m on the order of 0.5. Again, taking a five foot 
ius of curvature, the abscissa of Fig. 7 would be 
, and it may be seen that the heating rate is some- 
at over 100 BTU/ft?/sec even at low re-entry angles. 
would appear, then, that from the heating rate 
ndpoint, R CpA/m is going to have to be greater 
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Fig. 7. Design Parameter Diagram, Satellite Re-entry 
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than 2.5, and the larger the better. If it could be made 
as large as 200, this would eliminate the need for special 
heat protection devices. 

The peak deceleration starts at a little less than 8 g’s 
at shallow angles, and increases more or less linearly 
with angle to nearly 30 g’s at 10°. For manned flights, 
the maximum allowable deceleration probably lies 
somewhere between 10 and 12 g’s for the time dura- 
tions involved. This means that the steepest allowable 
angle would not be much over 3°. 

Impact uncertainty follows the opposite trend. 
Impact becomes more accurate as re-entry angle is 
increased. For the manned case, observe that 10 g’s 
peak deceleration corresponds to an impact uncer- 
tainty of 40 miles under the assumptions used. It should 
be emphasized that this 40 miles does not represent the 
total error. It represents only the error which results 
from uncertainties in re-entry angle, re-entry speed 
and atmospheric density uncertainties. There are other 
uncertainties to be considered, principally the point 
on the trajectory at which the vehicle crosses the re- 
entry altitude. This is determined by what happens 
outside the atmosphere, and is thus beyond the scope 
of study. This will probably be the largest single source 
of error. 

We see, then, two opposing tendencies: to get accu- 
racy the re-entry angle must be steep, to keep down the 
peak deceleration the re-entry angle must be shallow. 
Just how the compromise is to be made depends on 
the mission of the vehicle. 


SEC IV—Re-entry from Satellite Orbit. No 
Lift, Variable Drag Coefficient 

One method which has been proposed for alleviating 

the re-entry problem is that of varying the drag during 

the descent. It was indicated in the preceding section 

that so long as the drag was constant during the run, 
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the peak deceleration was the same, regardless of what 
the drag was. Heating rate and dynamic pressure did, 
of course, depend on CpA/m. If the drag is changed 
during the re-entry, however, this is no longer the case. 
It is possible to achieve some reduction in peak de- 
celeration. The purpose of this section is to show how 
much the deceleration may be reduced, and to investi- 
gate the effect on other parameters. 

Assume that the vehicle has a means of varying the 
drag. Let us define the maximum drag by CpA/m, 
and the actual drag by (CpA/m)e, where ¢ is a variable 
factor which lies on the range 0 < « S$ 1. The velocity 
equation for the flat, non-rotating earth with gravity 
neglected would be, then, 


dV CpA 
SSS ‘ 34 
tLe ( z ) : (34) 
Again using the exponential atmosphere, 
dV PO +72 Cp A —Kh 
oe a smn See ‘ 35 
dt 2 3 ( m \e ee) 


If we now apply the transformations as before: 


- ah waived AS oe et ten 
V=Vev, «=Kh (S25) (36) 


tKVegsngd =r 
then Eq 35 becomes 


Oa wate. (37) 
The optimum way to make use of variable drag would 
be to start the re-entry with the drag a maximum 
(e = 1). It would be maintained at this value (in order 
to keep the heating down) until the deceleration has 
reached the maximum desired value. Then ¢ is varied 
so as to keep the deceleration at this desired maximum. 
It is this portion of the process which must be con- 
sidered now. Generally speaking, if the deceleration is 
to be held constant as the altitude decreases, then, 


i (vee *) =0 = —ve+v—-+2 bgt) (38) 
dx dx dx 


It may easily be shown that in this approximate 
problem, 


—— = vee. (39) 


—=e-—22¢” (40) 
dx 


which may be combined with Eq 39 to give 


i v 
ies op yl ial (41) 
€ v 
This may be integrated directly to give 
Ine =  — 2Inv + Const. (42) 
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PEAUW DECELERATION , VARIABLE DRAG 
PEAK DECELERATION, CONSTANT DRAG 


Fia. 8. Effect of Drag Variation I 

It is now necessary to apply the initial conditions 
Recall that the above differential equations describ 
the portion of the trajectory where the deceleration i 
being held constant. The initial conditions for this are 
then, those of the point of the trajectory at which i 
was decided to begin holding the deceleration constant 
Assume this occurs as an x of x», and a v of vo, then 


et 7) 


() 4 


dv Gronks 
=~ = U = 00) = 
dx y Pe” =—eo 
— v 
Vo 
vd) = vee” ax. (44 


Integrating this equation, and applying the initia 
conditions, gives the results that 


— = 1+ 2¢°(2 — a). (45 
on 
Combining this with Eq 43 gives 


r—T 


e 
1 + 2e-*9(% — x) 


By differentiating this expression with respect to « ant 
setting the derivative equal to zero, the minimum valw 
of « may be obtained. It is 


evo 
: (1 ie ) (47 


€min > ais 


2e-*0 


(< 


(46 


This may easily be computed as a function of 2» 
From Table 1 it may be determined what the decelera 
tion is at each value of x). From a cross-plot of th 
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Fic. 9a. Normalized Deceleration vs Normalized Altitude 


two, the enin corresponding to each deceleration reduc- 
tion may be shown. This is given in Fig. 8, as the solid 
curve. 

_ Figure 9A is a plot of normalized deceleration as a 
‘unction of altitude, and Fig. 9B is a plot of € as a func- 
a of altitude. The deceleration reductions as taken 
rom curves like these have been plotted along with the 
analytical result in Fig. 8. Three cases were studied, 
2°, 3° and 5° re-entry angles. It may be seen that the 
© case agreed with the analytical result the best 
ecause, as mentioned earlier, the approximation is 
etter for this case. The other angles disagree some- 
hat more. 

It may be seen from Fig. 8 that it is a relatively easy 
atter to reduce the peak deceleration by 10, 20 or 
30 per cent. Any more than this begins to come hard, 
owever. In order to get a 50 per cent reduction, it 
ould take a rather ambitious drag varying device, 
mething more than 10 to 1. The peak heating rate 
s unaffected by this drag variation. This is because 
t least for all the cases studied here, the peak heating 
ad occurred before the drag started to vary. There is 
ome additional heating at lower altitudes, and the 
otal heat will be somewhat higher, but this is not really 
ppreciable. The impact prediction accuracy is not 
significantly affected. 

Various simplifications of this scheme suggest them- 
elves. The first is that of not increasing drag again 
fter it has passed its minimum. The deceleration would 
tart falling immediately after the minimum was 
assed, and there would be some additional heating at 
ower altitudes after the maximum had passed, but 
here would be no great difference, and the primary 
bjective of keeping the deceleration below the pre- 
eribed maximum would be attained. 

Another possibility is that of a sudden stepwise 
eduction in drag rather than the gradual, programmed 
eduction shown in Fig. 9. This was investigated and 
und to be incapable of reducing the peak decelera- 
ion. The deceleration dropped immediately, of course, 
vhen the drag was reduced, but later it built up again 
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Fic. 9b. Drag Variation vs Normalized Altitude 


to higher values, so that it made little or no reduction 
in the peak deceleration. 

In summary, then, it may be said that drag variation 
offers a possibility of modest reductions in peak de- 
celeration without significant penalty in any other 
area except for the additional complexity of the drag- 
varying device. If sizeable reductions in deceleration 
are desired, however, some other device will have to 
be found. 


SEC V—Lifting Re-entry from Satellite Orbits 


It has been proposed to use lift to alleviate various 
re-entry problems and, in fact, it shows real advantages. 
Assume that we apply a constant C,/Cp throughout 
the re-entry. Figures 10A, B and C show the effect of 
this lift on peak deceleration, peak heating rate and 
total heat input, respectively. Results are shown for 
each of two re-entry angles: 2 degrees and 5 degrees. 
The reduction in deceleration is approximately the 
factor 5. The improvement in peak heating rate is not 
so remarkable, being slightly less than the factor 2. It 
may be seen that the 5 degree re-entry case is less 
affected by the lift. This is because, with the steeper 
re-entry angle, the lift has less time to curve the flight 
path before peak heating and peak deceleration occur. 
It might be assumed that even better results could be 
obtained with shallower re-entry angle. 

So far, the use of lift appears to be quite attractive, 
but its effect on impact accuracy is somewhat more 
troublesome. Possibly Fig. 11 will give some idea of the 
nature of the difficulty. It may be seen that as the lift 
goes up, more and more ‘“‘skips” take place. Now these 
skips are somewhat difficult to control. Once the lift 
has been applied, and the vehicle starts up, it goes into 
a pure ballistic phase, and it is not possible to control 
it again until it comes back down into the denser 
atmosphere. Thus all points between the point it 
leaves the dense atmosphere and the point it returns 
are inaccessible to the missile as impact points. The 
length of these skips is rather critically dependent on 
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Fig. 10a. Effect of Lift on Peak De@ 
celeration. 
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the amount of lift applied. On all the trajectories, a 
point is finally reached where the missile has lost most 
of its initial speed, and it goes into a nearly vertical 
fall, at its terminal velocity. This is the point which 
must be controlled, because relatively little control is 
available thereafter. Let us say that h; is the altitude 
at which the vehicle loses all its forward speed (aside 
from that which might be supplied by lift), and call 
AR the amount by which the impact point may be 
changed by lift. It may be easily shown that 
iN tee Ore 
ics 

In no practical case will h; be greater than about 180,000 
feet. Assuming C,/Cp of unity would mean that the 
range could be changed by 180,000 feet in either direc- 
tion, or about 30 NMI, hardly enough to wipe out the 
effect of some erroneous skips earlier in the trajectory. 
In order to use lift to control the impact point, it 
appears necessary to use a rather elaborate tracking 
and computer system on the ground and to transmit 
commands to the missile from time to time during the 
re-entry. This is a complication which should be 
avoided if there is any other way to solve the heating 
and deceleration problems. 

In Fig. 11, the open circles indicate the point on each 
trajectory where maximum deceleration occurs, and 
the crossed circles indicate that point where maximum 
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Fig. 10b. Effect of Lift on Peak Heat- 
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Fic. 10c. Effect of Lift on Total Peak. 


heating occurs. Observe that on all the trajectorie 
the maximum heating occurs early, not long after the 
lift has begun to separate the paths. The maximum 
deceleration, however, occurs just before the final 
plunge into the dense atmosphere. This further il- 
lustrates the need to use the lift for deceleration re- 
duction until late in the trajectory, and only then is 
it free for use in impact control. 7 

A somewhat different situation exists for the steeper 
re-entry angle as is shown in Fig. 11. Here, both the 
peak heating and deceleration occur early in the tra- 
jectory, even before the first skip. The impact accuracy 
is much better at the steeper angle, but, on the othe 
hand, with the steeper angle, lift does not offer sé 
much improvement in peak deceleration and heating 
rate. i 

In summation, then, it appears that from the stand- 
points of deceleration and heating rate, lift is very 
attractive. It raises such problems in impact control, 
however, that it should be used only if simpler schemes 
prove unworkable. 


SEC VI—Re-entry from Lunar Orbits. No Lift. 
Constant Drag Coefficient 


As the next step in the study, it was decided te 
proceed to a higher-speed re-entry, consistent with @ 
circumlunar orbit. Only the concept of entering the 
atmosphere on the first pass was considered in any 
detail. The intent of the present study is to outline the 
problems of entering directly. That is to say that 
once the vehicle enters the atmosphere, it stays in 
until impact at the surface. 

Only one re-entry speed was studied in this case, 
that of 36,000 ft/sec. Figure 12 shows the character: 
istics of re-entry at this speed. Figure 12A is the peak 
deceleration as a function of re-entry angle. This curve 
has a minimum near 5°, and the curve does not g¢ 
below 4.8 degrees or so. The reason for these thing: 
may be gleaned from Fig. 13. This is a trajectory 
plot, very similar to that of Fig. 11. For the steeper 
re-entry angles, the trajectory proceeds more or les§ 
directly into the denser atmosphere—as might be 
expected. At shallower angles, however, the behaviol 
is different. Consider, for instance, the 5° trajectory 
The flight path angle decreases after re-entry starts 
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- Fie. 12a. Peak Deceleration vs Re-en- 
try Angle for Re-entry Speed of 36,000 
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until it finally reaches zero, and the vehicle is moving 
along the local horizontal. Because of the considerable 
drag existing at this altitude, it finally falls below the 
satellite speed and starts dropping into the denser 
regions. It is during this final plunge that the maxi- 
mum deceleration occurs, though the maximum heating 
occurred shortly after re-entry. Because of the high 
vehicle velocity, there is a strong tendency for the 
vehicle to ‘“‘graze’”’ the atmosphere, and leave again. 
This is exactly what happens if the re-entry angle is 
much less than 5°. The 4°50’ case does in fact leave 
the atmosphere again, as shown in Fig. 13. For refer- 
ance, the trajectory for a 5 degree re-entry angle is 
shown when the effect of atmosphere is deleted. The 
point of each trajectory where maximum deceleration 
xecurs is indicated by the open circle, and the poimt 
where maximum heating rate occurs is indicated by 
the crossed circle. 

Figure 12 shows also the values which would be 

predicted by the simplified analytical method of Sec. 
[IC. It may be seen that in all these cases, the analytical 
solution is in error by something like the factor 2. 
[his is probably due primarily to the fact that the 
arth cannot be considered flat for these trajectories, 
ind this was one of the assumptions in the analytical 
levelopment. The agreement between simulator and 
malytical results was relatively good for satellite 
e-entry but it is quite poor for lunar ones. 
- Figure 12B shows the peak heating rate. There is no 
ninimum on this curve because peak heating always 
ecurs early in the trajectory before curvature of the 
arth or anything else has had much effect. Figure 12C 
hows the total heat input, and it is what would be 
xpected from Fig. 12B. 

Just as in the satellite case, the next consideration is 
hat of uncertainty in the impact point. Range versus 
e-entry angle for 36,000 ft/sec re-entry speed is shown 
n Fig. 14. Figure 15 shows the sensitivities, and I'ig. 
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Fig. 12c. Total Heat Parameter vs 
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16 the errors which would result, using the same un- 
certainties as in the satellite case: 30 per cent density 
uncertainty, 0.1° angle error, 20 fps velocity error. 
The impact uncertainties run considerably higher 
here than for satellites. 

All these considerations are summarized in Fig. 17, 
which is analogous to Fig. 7, and is used in just the 
same way. In fact it is quite interesting to compare 
the two. First of all, it may be seen that the lunar 
heating curves are shifted to the right by about the 
factor 5. This means that to get the same heating, for 
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Fig. 13. Lunar Re-entry Trajectory Plots 


instance, the lunar RCpA/m must be about five times 
greater than the one for satellite re-entry. Putting it 
the other way around, if the same RCpA/m is used in 
both cases, the lunar heating rate will be about /5 = 
2.2 times greater. Of course this depends on the re- 
entry angle used in the two cases, but it is indicative 
of the difference. It is interesting that the peak decelera- 
tion for the lunar case may actually be somewhat less 
than the best obtainable in satellite re-entry. Possibly 
the most significant difference between the two cases 
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Fig. 14. Range vs Re-entry Angle for 
Re-entry Speed of 36,000 ft/sec. 


is in impact point uncertainty. The lunar ones are much 
higher. Assume for instance, that we wish to limit the 
decelerations to 10 g’s in both cases. The satellite un- 
certainty corresponding to this is about 40 NMI, 
while for lunar re-entry, it is more like 100 NMI. 
In the latter case the uncertainties are such as to cast 
doubt on the feasibility of a non-lifting re-entry. 


SEC VII—Re-entry from Lunar Orbits. Low 
Lift, Variable Drag Coefficient 


It is possible to use variable drag to decrease the 
peak decelerations in lunar re-entries, just as in the 
satellite cases. The same program may be followed, and 
with much the same result. The reduction in peak 
deceleration as a function of ¢, the drag-variability 
factor falls exactly on top of the 2° curve of Fig. 8. 
It does not affect peak heating rates, increases total 
heating slightly, and does not affect impact accuracy. 

Lift, in this case is not nearly so effective as in satel- 
lite re-entries. It was found that if lift is applied early 
in the trajectory, in the region of maximum heating 
and maximum deceleration, then it is not possible to 
apply more than C,/C of more than 0.1 or 0.2 without 
causing the vehicle to leave the atmosphere altogether. 
Applying lift later in the trajectory, of course does not 
decrease the heating or deceleration, though it provides 
a means of changing the impact point by a considerable 
amount. Again, the question of accuracy in controlling 
C,/Cp comes into play, however. The accuracy re- 
quirements are about the same as those given earlier 
for the satellite case. 


SEC VIti—Results and Conclusions 


First of all, it is of interest to present a comparison 
of the lunar and satellite re-entry problems with each 
other, and with the more familiar ballistic missile 
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re-entry problem. Figure 18 gives some aspects of thi 
comparison. In all three figures, the re-entry angle i 
plotted as the abscissa and a different scale is use 
for the ballistic missile parameters from that used fo 
satellite and lunar curves. Figure 18A shows the peal 
deceleration comparison. For larger angles, more thai 
100 g’s of deceleration is possible. Satellite and luna 
decelerations fall in the same range, from 8 to 20 g’s 
and both are well below the corresponding figures fo 
ballistic missiles. From this standpoint, the satellit 
and lunar re-entries are easier, and this is principalk 
because of the much shallower re-entry angles. Figur 
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Fig. 17. Design Parameter Diagram, Lunar Re-entry 


18B shows the peak heating rate comparison. Her 
both lunar and satellite values lie above those for th 
IRBM and below those for the ICBM, and the luna 
values run considerably higher than the satellite one: 
It should be kept in mind that these comparisons at 
sume the same values of R and Cy)A/m for all the classe 
of vehicles. This may or may not be the case, but th 
assumption seems to provide the only meaningft 
basis for comparison. 
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Figure 18C shows the comparison of total heat input, 
id it is here that differences are most marked, and 
1e of the major difficulties in satellite and lunar 
-entry is indicated. The satellite figures run three 
mes higher than the ballistic missile ones, and lunar 
-entry some five times higher. As explained earlier, 
is is total aerodynamic heat input, and says nothing 
yout how much is re-radiated. In the steeper angle 
ses, little is re-radiated, though for shallower angles, 
-radiation may be an appreciable amount of the 
hole. In fact if CpA/m is high enough, all the in- 
ent heat may be re-radiated. 

In regard to impact accuracy, of course, ballistic 
ssiles are far superior to the other vehicles. The 
in reason is the much steeper re-entry angles used. 


Satellite Re-entry Discussion 


The simplest system for re-entry is a non-lifting, 
istant-drag vehicle, which follows a simple ballistic 
th. Let us consider the three areas of major concern: 
celeration, heating and impact error. The peak decele- 
ion can be held to 8 g’s with large impact errors, to 
g’s with more reasonable impact errors. Impact 
curacy may be disposed of by saying that it appears 
ssible to achieve accuracy sufficient that proper 
und installations, suitably alerted, could find the 
hicle. It does not seem that it will ever be possible 
th this simple system to deliver the vehicle with 
uracies sufficient for warhead delivery. If this is 
be done, some additional control must be provided. 
re will be said about this presently. 

eating is probably the most serious problem. In 
er to avoid ablation or heat-sinks, it would appear 
essary to keep RCyA/m above about 150. This 
uld mean, for instance a 5 ft radius of curvature and 
',A/m of 30. If excessive ablation rates are to be 
ided, RC,»A/m must be kept above about 4. With 
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a 5 ft radius of curvature, this would mean a CpA/m 
of 0.8. If a heat-sink is to be used, it would have to be 
large compared to what is used in ballistic missiles, 
because of the much larger total heat input. BTU/ft? 
of 5000 to 15,000 would be possible for reasonable 
designs. Study of the heat protection question itself is 
beyond the scope of the present work. The only intent 
here is to show the nature of the problem. 

Let us now consider possible improvements in ac- 
curacy. In a satellite recovery, the main point at 
issue is the predictability of the impact point. Since 
the satellite is presumably subject to commands at 
any point of the orbit, the point at which the recovery 
is initiated may be easily controlled by the designer. 
The requirement, then, is to predict how far from this 
initiation point impact will take place. This is in 
contrast with the lunar case, where there is no flexi- 
bility whatever in the point at which recovery is 
initiated, and it may be desirable to vary the length 
of the re-entry trajectory over a wide range. To return 
to the satellite case, since predictability is the point, 
and not wide-range control, any mechanism which can 
eliminate the 100 or 200 NMT of error which would 
normally exist, will be adequate. From Fig. 6 it may be 
seen that the principal source of error is uncertainty 
in the atmospheric density. Next comes re-entry angle 
errors. If really precise impact point control is desired, 
some way will have to be found for overcoming both 
these error sources. It would appear that drag variation 
might do it. The range could be changed by 200 miles 
or so by changing the drag during descent. This would 
have to be done with some care to avoid heating and 
deceleration problems. This area requires further work. 
In principle, lift could be used to do the same job, but 
it appears, at least on the surface, that this would be 
difficult to accomplish, because lift is so difficult to 
measure. This too, must be studied further. 


The Journal of the Astronautical Sciences 19 


B. Lunar Re-entry Discussion 


In this case, deceleration does not appear to be a 
major problem. It is possible to get decelerations even 
lower than those for satellites, though at a severe 
penalty in impact accuracy. 

Heat protection is critical, and now so is impact 
accuracy. In this case there is little that can be done 
to improve the heating situation except to abandon 
the idea of re-entry on the first pass. Lift or drag 
variations seem to offer little improvement. To dissipate 
all heat by radiation, it will be necessary to keep 
RCyA/m above 600 or so. To avoid excessive ablation 
rates, it would have to be above 25. Even assuming a 
10 ft radius of curvature, this would mean CpA/m’s 
of 60 and 2.5 respectively. The first would appear to 
be out of the question. 

As for impact accuracy, the simple system is marginal 
at best, and this is leaving aside the not inconsiderable 
problem of guiding the missile to satisfactory re-entry 
conditions. From then until impact, the errors de- 
veloped are serious, and it seems likely that a simple 
ballistic system will be unfeasible. This is all from 
the standpoint of predictability alone. It seems quite 
desirable, and in fact almost a necessity to be able to 
vary the re-entry trajectory over rather wide limits. 
If this is not done, the impact point will be determined 
by the orientation the earth happens to have when 
the vehicle approaches. While this could be predicted, 
and accounted for at launch, some flexibility is desirable. 
There will be various other factors to be considered in 
selecting a desirable launch time, and it may not be 
possible to satisfy all requirements simultaneously. 
There are two obvious ways of supplying this flexibility. 
One is the use of “braking ellipses” to slow the vehicle 
down in several passes into the atmosphere, until a 
more or less circular orbit is achieved, from which 
final recovery may be initiated. The other method 
would be to enter on the first pass, but to supply lift 
or other control after the vehicle had fallen below 
satellite speed. Further study will be required, but 
based on present evidence, it is felt that the second 
alternative is the more preferable. Braking ellipses 
aggravate the guidance problem and lift inside the 
atmosphere aggravates the total heat problem though 
it would not increase the peak heating rate. Neither 
of these problems have yet been solved, so it may be 
idle to speculate on which is the more impossible. 
High temperature sublimation and ablation materials 
have appeared on the horizon, whereas the same cannot 
be said of lunar guidance equipment. 

Of course it will be necessary to supply guidance to 
reach the correct re-entry point in either case. For 
instance, the listed tolerance of 0.1° in the re-entry 
angle corresponds to a difference of only 6,800 feet in 
the perigee which the orbit would have in the absence 
of any atmosphere. In other words, the guidance 
system would be required to control the hypothetical 
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perigee or distance of closest approach to this ai 
curacy. The perigee would occur, incidentally belo 
200,000 feet, depending on the re-entry angle usec 
For braking ellipses, it would appear that this accurae 
would have to be supplied, not once, but several time 
and uncertainties in atmospheric density would f 
quite serious. 


Nomenclature 


D Drag 

K Constant, Equation (2) 

K Reciprocal of Body Nose Radius 
L Lift 

Lz Lewis Number 

Q, Aerodynamic Heat Input 

Q. Aerodynamic Heating Rate 

R Body Nose Radius 

R, Mean Earth Radius 


Ag 
: 
+ 
: 
4 


g Gravitational Acceleration at Earth’s Surfac 
h Altitude above the Earth ; 
h Gas Enthalpy 
hp Average Dissociation Energy per Unit of Mas 
7 Joule’s Constant 
k Density Ratio Across the Shock 
q Dynamic Pressure 
ous Drag Parameter 
m 
Bus Lift Parameter 
m 
U Velocity 
2h eS oR) Non-dimensional Stagnation Poin 
OG \Ga ys Gradient 
¥ Flight Path Angle 
6 Negative Flight Path Angle 
ro) Range from Re-entry 
p Air Density 
m Viscosity Coefficient 
) Shock Layer Thickness 
Subscripts 
s Stagnation Point 
Ww Body Wall 
0 Free-Stream 
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Technical Notes 


Dual Burning Propulsion Sys- 
tems for Satellite Stages 
bei Marti 


.bstract 


_ -A satellite vehicle system whose last stage can be shut 

down and restarted is shown by basic energy relation- 
-ships to have a decided performance advantage over 
a vehicle system which is identical except for having 
only the conventional one-shot capability. For an as- 
sumed final stage weight and boost ascent, curves of 
payload vs. orbit altitude are derived and compared for 
the two final stage thrust systems. 


troduction 


‘Except for effects of gravity and drag, a given impulse 
rected along the velocity vector of a mass will result in the 
ume velocity change of the mass for all values of original 
plocity of the mass. But since kinetic energy varies as the 
yuare of the velocity magnitude, the energy change for a 
ven impulse increases with increasing values of original 
locity. There is, therefore, an advantage, from the stand- 
bint of minimizing rocket propellant required, in supplying 
much of the total impulse as possible when the velocity 
the rocket is high; i.e., when the energy of the rocket is 
rgely in kinetic form. Ideally then, all of the energy re- 
ired for a particular satellite orbit should be supplied just 
bove the atmospheric drag regions, whereupon the vehicle 
oves on a frictionless track into a horizontal trajectory at 
e proper altitude. Since these frictionless tracks are not 
hder development, the next best technique must be em- 
oyed. In this technique, as little impulse as possible is left 
be supplied under conditions of reduced velocity and high 
itude and is accomplished by having the begin-coast 
itude as low as practical and the begin-coast velocity as 
gh as possible consistent with the final altitude desired. 
he limiting value for this velocity is that for perigee of an 
icent ellipse whose apogee altitude is essentially that of the 
»sired final orbit. 

wal Burning Concept 

For good staging ratios, however, a stage break probably 
ill not be indicated at perigee, in which case only the apogee 
blocity increment would be supplied by the final stage. 
|* Staff Scientist, Spacecraft Dept., Lockheed Missiles & 
pace Division. 
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ASSUMPTIONS 
BEGIN-COAST ALT.=100(NM) 


VELOCITY AT 1OO(NM) DUE TO 
BOOSTER = 18,000 FT/SEC 


FOR DUAL-BURN, COAST RANGE 


SATELLITE STAGE GROSS WT 
=20,000 LB 


lap * 350 SEC 
FIXED INERT WT = 2000 LB FOR 
SINGLE-BURN 6 2100 LB FOR 
DUAL-BURN 
PROPELLANT SYSTEM WT « 
4% OF IMPULSE PROPEL- 
LANT WT 


PAYLOAD WEIGHT (LB) 
ro) 
° 
to} 
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SINGLE-BURN 


ORBIT ALTITUDE (NM) 


Fie. 1. Payload vs. Orbit Altitude 


For most Earth satellite orbit altitudes and vehicle con- 
figurations, high staging efficiency will call for the final stage 
to assume a greater part of the overall velocity requirement 
(frequently referred to as characteristic velocity) than the 
velocity increment required at apogee. These ascent and 
performance criteria led to consideration of two periods of 
thrust for the final orbiting stage; 1.e., “dual burning’’, in 
order to arrive on the coasting ellipse near perigee while an 
efficient stage distribution of the total velocity requirement 
is maintained. 

The more conventional single-burn trajectories and dual- 
burn trajectories may be summarized for comparison as 
follows: 

(1) Single-Burn Ascent Trajectory—At burnout of the 
boost stage, or stages, the flight path angle is such 
that for the speed and altitude attained, the coast 
phase which follows will carry the vehicle to apogee 
which is essentially equal to the final orbit altitude. 
Just prior to apogee the final stage propulsion system 
is operated to provide orbit injection. 

(2) Dual-Burn Ascent Trajectory—At burnout of the boost 
stage, or stages, the flight path angle and altitude 
are, in general, lower than for the single-burn case. 
Immediately after booster separation, the final stage 
propulsion system is operated long enough to in- 
crease the speed essentially to perigee speed for an 
ascent ellipse whose apogee is at the altitude of the 
desired orbit. After a long coast (through a range 
angle of up to 180 degrees) the propulsion system 1s 
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restarted and operated long enough to bring the 
speed up to orbital speed. 

As mentioned earlier, total velocity requirements are 
lowest for a coast range of 180 degrees (coast begins at 
perigee) because this permits the highest begin-coast speed. 
However, the effect on the total velocity requirement of 
reducing coast range from this value is rather slight until 
coast range is substantially reduced. Therefore, the exact 
coast range for a particular mission should be selected after 
including other considerations such as guidance component 
errors for various coast times, sensitivity ranges of radiation 
sensors if used in the ascent guidance system, and propellant 
boil-off during coast. For the performance comparison of 
single-burning and dual-burning which is presented herein, 
a coast range of 180 degrees was assumed for the dual-burning 
case. 

For both cases energy is conserved during the coasting 
ellipse, so that the vehicle velocity at the end of the coast 
period may be found from the kinetic and potential relation- 
ship as follows: 


CE i en Bg (1) 


where: V; = velocity at beginning of coast 
r, = length of radius vector at beginning of coast 
GM =the product of the universal gravitational 
constant into Earth’s mass, taken as 1.4077 
< 10"? ft? sec 
V2 = velocity at end of coast 
ro = length of radius vector at end of coast 


Vo = / Vie eee (2) 
T2 Tr 


It was assumed for the example that the satellite stage gross 
weight and the booster vehicles are identical for both types 
of trajectories, and that the mission calls for a circular orbit 
at 1000 nautical miles altitude. 

Booster burnout for the single-burn vehicle was assumed 
to occur at 100 nautical miles altitude at a velocity of 18,000 
ft/sec and for the dual-burn vehicle, at a somewhat lower 
altitude (due to the flatter trajectory) but with the same total 
energy. 


or 


Single-Burn Performance 


After booster burnout, the single-burn vehicle coasts on a 
constant energy ellipse to apogee at 1000 nautical miles. 
The begin-coast conditions coincide with the booster burnout 
conditions. The velocity at the end of the coast period is 
found from Kq. 2 as: 


V2 = 7650 ft/sec. 
The required orbital velocity at 1000 nautical miles is 
22,850 ft/sec. The velocity differential, or Av, that must be 


supplied by the vehicle then is 22,850 — 7650 = 15,200 
ft/sec. 


Dual-Burn Performance 


In the case of dual burning, the boost phase is immediately 
followed by the first of two separate thrust phases. The 
desired velocity (perigee velocity) at the end of this first 
burning phase for a coast range of 180 degrees may be found 


from: 
Vp = 4/ 50M) Senses (3) 
rp(rp + 14) 
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where: Vp = velocity at perigee 
r4 = apogee radius = Earth’s radius + 1000 n.m. 

rp = perigee radius 

The value of rp depends upon the individual vehicle an 
trajectory characteristics. For this example assume that 
altitude at the end of the first burning phase, and the beg 1 
ning of the coast phase (perigee), is again 100 nautical miles 
Since the total energy at booster burnout was assumed to b 
equivalent to that for single burning, the velocity incremen 
for the first burning phase can be taken as the differene 
between the perigee velocity and 18,000 ft/sec. 
This is a rather conservative approach for two reasons 
viz.: (1) Since the boost trajectory will be somewhat flatte 
for dual burning, gravity losses will be lower and burnou 
energy will be higher in spite of slightly increased drag losses 
(2) the first burning velocity increment is actually accumu 
lated over the altitude change from booster burnout to begii 
coast which (because the impulse is furnished at a mor 
favorable kinetic-potential energy distribution) is mot 
efficient than the assumed instantaneous burning at begit 
coast altitude. 
For the assumed 100 n.m. perigee, perigee velocity, fror 
Eq. 3 is: 


Vp = 26,990 ft/sec. 

The velocity increment to be supplied during first burn i 
then: : 
' 


AV, = Vp — 18,000 = 8990 ft/sec. 


The velocity after coasting to 1000 n.m. may be fount 
(by conservation of energy) from Eq. 2, or, by conservatiol 
of angular momentum: 


Va = Volrp/r a] 
Va 21,520 ft/sec. 


(4 


I 


The velocity increment required at apogee during thi 
second burning period is then: 


AV» = 22,850 — V4 = 1330 ft/sec. 


The total satellite stage velocity increment required fo 
dual burning is: 


AV, a AV = 10,320 ft/sec. 


Payload Comparison for Single-burning and Dual-burn 
ing Stages 


The above velocity requirement of 10,320 ft/sec. is com 
pared with 15,200 ft/sec. required for single burning. Th 
difference in final stage velocity requirements can be con 
verted into payload differential for a given final stage gros 
weight and known characteristics of the propulsion systen 
and other subsystems. This velocity-to-payload conversiol 
was computed for the boost ascent assumed herein for severa 
orbit altitudes based upon a satellite stage gross weight @ 
20,000 Ib. and specific impulse of 350 seconds. Payloas 
weight was assumed to be related to burnout weight by th 
following expression: 


Wet = Wao — 2000 — 0.04Wp — X @ 


where: Wp, = weight of payload 
Wo = weight at burnout = 20,000e~‘4¥ /35o” 
Wp = weight of impulse propellant required 
X = zero for single-burning and 100 lb. for dual 
burning to cover second start requirement 


ll 


The results are compared in Fig. 1. The boost trajectory 
d payload expression, although greatly simplified, serve 
show that substantial payload increases are available. 
e abrupt decrease in payload for the single-burn curve is 
e to the rapid increase (to vertical) for the begin-coast 
ht path angle to achieve the indicated altitude at apogee. 


The maximum altitude obtained for the above assumptions 
with single burning was found to be 1265 nm. For this 
altitude, since the ascent must be vertical, apogee velocity 
is zero and the velocity increment required is the entire 
orbit velocity. 


Anisotropy of Escape Velocity from the Moon, the Lunar 
Atmosphere and the Origin of Craters 
Louis Gold* 


pstract 


A simple estimate for the degree of anisotropy to be 
expected in the lunar escape velocity owing to the earth’s 
gravitational field provides a basis for assessing the role 
of directed streaming of the moon’s atmosphere—the 
lunar wind—as it might have contributed to the genesis 
of craters. The constant aspect of the moon (relative to 
the earth) might have provided the essential orienting 
effect in the formation of the ordered appearance of the 
lunar craters during its primordial evolution. The pos- 
sible connection of the streaming ionized residual at- 
mosphere of the moon with a predicted weak but 
decidedly anisotropic lunar magnetic field is pointed out. 
Thus the ideal launching positions for return flights to 
earth may have to allow for dangerous radiation belts 
encircling the moon as well as minimal propulsion re- 
quirements, 


scussion 


It was pointed out some time ago [1] that the strong ani- 
rropy of the escape velocity from the surface of the moon 
ght have a bearing on the formation and general ap- 
arance of the lunar craters. This brief note is intended to 
‘borate somewhat on such an hypothesis and its con- 
juences. The imminence of lunar landings and exploration 
ikes it of practical import to report on the following theo- 
ical considerations. Recent accounts of volcanic activity 
the moon have also stimulated current interest in the 
ystery of the lunar craters [2]. 
While the restricted three body problem as solved in the 
m of the Jacobi integral does imply an anisotropy of the 
‘ape velocity from the moon, this has not been generally 
preciated. Moreover, no actual calculation of this effect has 
parently been made. Since the Jacobi solution fails to 
corporate many other possible contributions to the lunar 
‘ape such as the asphericity of the earth and moon, solar 
rturbation, etc., the approximation derived earlier [3] 


igeeg M, M, M2 M 
a — v2) = G ( iz +3) ae G Got ) 


yy be employed for this purpose. Here at time t = 0 the 
tial velocity v = v and M,, M, denote the respective 
ses of the moon and earth; d is the distance between the 
ters of MW, and M2 with 7» the lunar radius, etc. The sign 
propriately must be chosen according to whether flights 
ginate on the near or far side. 

Escape for the latter case corresponds tov = Oasr— », 
ereby 


* Radiation, Inc., Research Division, Orlando, Florida. 
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in which gp represents the conventional gravitational constant 
referred to the moon’s surface. On the other hand, for the 
near side, the escape velocity Vz associates with r = rg = 
d(1 + M»/M,)—, whence 


M, M, M, i) 
Veins. = 24 (— — ee \ Se : 
(V 2) vs (= d= -) worl M, d+ =) (2b) 


If the true neutral point is introduced 


pL Cyl play et (Ba) 
M, 


a somewhat altered expression for (2b) arises 
(Vz2)w.s. = 2G a ae ) 
B)N.S. = 4 i 2p Rh 


1 1 
M, a 
at (, — ro +) 


) 


(3b) 


As only an approximate estimate is involved here anyway, 
the similarity of form in (2a) and (2b) permits easy insight 
as to how the ratio (Vz)n.s./(Vz)r.s. depends upon the 
masses and distances in the domain of interest. 
Numerically then the simple calculation 
Ms, To 


--—— & 0.34 
M d+ 


leads to the ratio 


(Views. és (; — 0.34\? = 0.70 
(J Seavey 1 + 0.34 


which demonstrates the significant role played by the earth 
in influencing flights leaving from the moon. But another 
factor which has to do with ideal launch from our satellite 
concerns the possibility of lunar radiation belts, an aspect to 
be entered into later. 

The severe anisotropy of the escape velocity which has now 
been characterized may have exerted an important influence 
on the origin of the lunar crater pattern. It may be conjec- 
tured that while the moon’s surface was in a semi-plastic 
state during its evolutionary phase with hot gas issuing from 
its interior, the severe streaming of the atmosphere produced 
huge whirls which were eventually frozen in upon final solidi- 
fication. The oriented streaming, maintained by the constant 
aspect of the moon presented to the earth, would tend to 
induce some sort of ordered appearance of the resultant 
craters: the meteor impact hypothesis hardly accounts for 
the general appearance. The craters of Tycho and Copernicus 
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are excluded from this picture as they may have been formed 
by meteor bombardment as suggested by the star-like rays 
surrounding them which have been ascribed to the settlement 
of meteoric dust in the moon’s gravitational field [4]. 

The primordial streaming must likely have been typical of 
mach flow in the hypersonic regime. While the present reced- 
ing position of the moon may considerably moderate the 
streaming action, it is possible that the greatly reduced 
density of the lunar atmosphere may yet be conducive for a 
high degree of ionization. Thus, apart from the question of 
concentrating the craters at the periphery of the lunar dise 
which corresponds to the observed situation, there is then the 
implication of a highly ionized residual lunar atmosphere. 
Failure to detect the moon’s atmosphere by close passage of 
radiation from a star may in some way be due to the pecu- 
liarities of such a streaming, ionized sheath. 

At this juncture, there arises the suggestion of highly aniso- 
tropic magnetic fields surrounding the moon; a Junar wind of 
the character described above has its evident counterpart of 
solar winds [5]. The prediction is made that the moon will be 


found to possess a weak but strongly oriented magnetic fiele 
The question of charge particle trapping naturally follow 
and so the possibility of lunar radiation belts may need to | 
appraised, particularly with regard to the contemplate 
manned landings. More detailed theoretical developmen 
along these lines are in progress. 
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Lists of figures, tables of contents, and distribution 
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Titles 


The title should be brief, but express adequately the 
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should be written as he prefers; all titles and degrees or 
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duction, covering contents of the paper. It should not 
exceed 200 words. 
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__ The paper can be divided into principal sections as 
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Mathematical Work 


As far as possible, formulas should be typewritten. 
Greek letters and other symbols not available on the 
typewriter should be carefully inserted in ink. Each 
such symbol should be identified unambiguously the 
first time it appears. The distinction between capital 
and lower-case letters should be clearly shown. Avoid 
confusion between zero (0) and the letter O; between 
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The level of subscripts, exponents, subscripts to sub- 
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Format of Technical Papers 
for AAS Journal 


Complicated exponents and subscripts should be 
avoided when possible to represent by a special symbol. 

Fractions in the body of the text and fractions occur- 
ring in the numerators or denominators of fractions 
should be written with the solidus. Thus 


cos (1x/2b) 
cos (ra/2b) 


is the preferred usage. 

The intended grouping of handwritten formulas can 
be made clear by slight variations in spacing, but this 
procedure is not acceptable in printed formulas. To 
avoid misunderstanding, the order of symbols should 
therefore be carefully considered. Thus 


: (a + bx) cos t is preferable to cos t (a + bx) 


In handwritten formulas the size of parentheses, 
brackets and braces can vary more widely than in 
print. Particular attention should therefore be paid to 
the proper use of parentheses, brackets, and braces 
(which should be used in this order). Thus 


{[a + (b + cx)"] cos ky}? 


is required rather than ((a + (b + cx)") cos ky)?. 
Equations are numbered and referred to in text as 
(15). 
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